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We study the Josephson effect in s±/I/s± junctions made by two bands reversed sign s-wave 
(s±) superconductive materials. We derive an equation providing the bound Andreev energy states 
parameterized by the band ratio a, a parameter accounting for the weight of the second band 
with respect to the first one at the interface. For selected values of the band ratio and tunnel 
barrier amplitude, we predict various features of the Josephson current, among which a possible 
high temperature ty state of the junction (a doubly degenerate junction ground state) and a ty — »• 
crossover with decreasing temperature. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.70.Dd 



The importance of the Josephson effect as a tool to 
probe the properties of new discovered superconductive 
materials can hardly be underestimated. This has been 
true in the past in investigating the d-wave cuprate ma- 
terials and is nowadays the case for the multi-gap iron 
based superconductors [1] whose complex behavior, be- 
yond the BCS theory, is presently a focus in condensed 
matter theory 0] and material science Q • 

There is a rather solid indication now, supported by 
experiments, that the pair potential symmetry in these 
compounds is s-wave with sign-reversing order parameter 
(s^). Several theoretical models have already discussed 
the experimental consequences of an extended s-wave 
(s ± - wave) order parameter symmetry on the Andreev 
conductance of an NS interface and on the Josephson 
effect in the iron-based superconductor junctions |5l-[l9|. 
In particular the Josephson effect has been studied in hy- 
brid devices i.e. sls± junctions, as summarized by Seidel 
in his review [20] and reference therein. 

In this Letter we study an all-pnictide symmetric 
s±Is± Josephson junction and we discuss a number non- 
trivial physical consequences on the Josephson effect due 
due to the presence of a second conduction band. We 
show that, depending on the band ratio parameter a, 
which accounts for the weight the second band, a 7r-state 
[2lj can develop for a wide range of junction transparen- 
cies and temperatures. This 7r-state can persist in the full 
range of temperatures or can undergo a tt — >• (inverse 
— » 7r) crossover as the temperature decreases, depend- 
ing on the band ratio parameter. This crossover is analo- 
gous to the — » 7r crossover [22, 23], found in mesoscopic 



d-wave superconductor Josephson junction (24[. 

We adopt the simplest model of a Josephson junc- 
tion that shows the essential features of the Josephson 
effect in the presence of two gaps, namely we consider 
a superconductor(S)-insulator(I)-superconductor(S) con- 
tact. The iron based junction is modeled by considering a 
one-dimensional conductor, whose left (x < 0) and right 
(x > 0) halves are both two band metals (two differ- 
ent states at the Fermi level, one with the wave vector 
p and the other with q). We assume that the motion of 
quasiparticles is described by the Bogoliubov de Gennes 
(BdG) equation [25| and that the order parameter has al- 
ready been obtained self-consistently from the gap equa- 
tion. Therefore we choose a one dimensional model for 
the gap, such that the left and right two band supercon- 
ductors have pair potentials given by 



Aj(x) = AjjwOi-x) + A j e i ^^ ) 9(x), j = 1,2 (1) 



The normal region, where Aj = 0, has an infinitesi- 
mal width and we also introduce a scattering potential 
U(x) = Uq5(x). The possibility of nodes in the gap func- 
tion is not considered. In the case of the two-gap model 
with unequal s-wave gaps, we write a wave function of the 
same t ype introduced in the Blonder Tinkham Klapwijk 
model [26] as solution of the BdG equations and treat the 
presence of the second band through the introduction of 
Bloch wave functions [8[. The bound state (B) eigenfunc- 
tion with energy \E\ < Ai (we assume Ai < A2) can be 
written as 
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where ip is a global phase difference between the two 
superconductive regions and </?i, </?2 are the phases of the 
gaps Ai, A2 in both p and g bands respectively. In the 
case of s± gap model, (p 2 — (fi = n - I n the considered en- 
ergy range the wave function has to decay exponentially 
for \x\ — > 00. The coefficients a B , b Bl c Bl d B are the prob- 
ability amplitudes transmission with branch crossing or 
without branch crossing. Essential above is the intro- 
duction of ao, a mixing coefficient defining the ratio of 
probability amplitudes for a quasiparticle to be transmit- 
ted to the first, (p), or second, (g), band [8]; the functions 
0's are the Bloch waves in the two-band superconductor; 
p and q are the Fermi vectors for the two bands corre- 
sponding to the same energy E [8[: 



M x ) = Yl Cg + exp ^ A + G ^ 



(3) 
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where A = p e ,Qe,Ph,Qh, and G represent any reciprocal 
lattice vector, m, v\ and u 2l v 2 are the Bogoliubov co- 
efficients for the first and second band, respectively 
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The global wave function \I> must satisfy the following 
boundary conditions at the interfaces x = 



^(0") = ^ 5k (0 + ) 



'S R m+\ M^S L m-\ - 2 y ^ g fl(Q+) (5) 



where primes denote derivative with respect to x. 
Coupling between the bands is implicit in the boundary 



condition requirements. In fact, the assumed wave func- 
tion \I/ is of the form \I/ = ^1 + ao^/ 2 where, separately, 
^1 and ^2 solve BdG equations for excitations of energy 
E for the gap Ai and A2, respectively,. The boundary 
conditions, Eqs. (j5j), are requested to be a constrain for 
the whole function ^. 

The above matching procedure, with the require- 
ment of non-triviality of the solution for the coefficients 
a Bl b Bl c Bl d B provides for the s±/I/s± junction, the 
spectral equation 



(-1 + 2yr^£Vr 2 -E 2 a 2 - r 2 a 4 ) (2Z 2 + 1 
2E 2 [Z 2 (l + a 4 ) + (l-a 2 + a 4 )] + 
"-1 + a 2 ( 2E 2 + 2y/l- E 2 ^r 2 - E 2 - r 2 a 2 J 
cos((p) = 



x 

(6) 



where we have introduced the gap ratio r = A2/A1, the 
band ratio parameter a = (0)/<j> p (Q) and the bar- 
rier strength Z = Uo/hv. The energy is given in units 
of Ai(T), and the two gaps will be assumed to obey the 
same BCS-like temperature law. Boundary conditions 
on the wave function derivatives are usually discussed in 
terms of Fermi velocities in the case of plane waves. For 
Bloch waves we have to introduce "interface velocities" 



m0 A (O) 



(7) 



In deriving the spectral equation we have assumed, for 
the sake of simplicity, equal band interface velocities v\ = 
v. 

An analysis of the spectral equation shows that for 
a 7^ there are, in general, four energy levels, E\ = 
d=ei(y?, a) and E 2 = ±62(^,0^) (see Fig. [T|. In a single- 
band case (a = 0) Eq [6] provides the well known result 
for a conventional SIS junction [27], HH, namely E\ = 

±e 1 (<p,0) = ±Ai(T) [1 - Dsin 2 ((^/2)] 1/2 , where D = 
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1/(1 + Z 2 ) is the transmission probability through the 
(5-function barrier, i.e. the junction transparency. 

As a increases, the energy levels E 2 = ±62 (y?, a) start 
to branch off the levels E\ = d=ei(y>, a), while these latter 
gradually approach the zero energy state. As an exam- 
ple, Fig. [1] (a) and (b) show the modifications of the 
Andreev levels for increasing values of the band ratio, 
for two values of the barrier parameter, Z = 0.7 and 
Z = 2 and for the band gap ratio r = 2, respectively. 
The condition for the existence of a zero energy level can 
be easily obtained from the spectral equation. It is given 
by (1 - ra 2 ) 2 (l + 2Z 2 + cos(y?)) = 0. Accordingly, a zero 
energy state E\ = ±ei(<p, a c ) = 0, is obtained for the 
critical value of the band ratio a c = 1/y/r, for any Z 
value. The energy level E2 = ±e 2 ((p,a c ) corresponding 
to the critical value of the band ratio is given by [29| 

E 2 ±e 2 ((p,a c ) 

1 2Dr 2 (2 - D + D cos cp) sin 2 (cp/2) 
Y 1 + r(2 - 2D + r) + 2Dr cos cp ^ 

and it is shown as a blue line in Fig. [I] for the indicated 
values of D and r. For low transparencies (D <C 1), 
Eq. [8] reduces to E 2 = ±e 2 ((p,a c ) = ±A 1 (T)2r/(l + 
r)y/~D sin ((/?/ 2), a result closely resembling the midgap 
states of d-wave superconductor Josephson junction [30|. 
For a > a c , there are no surface bound states with real 
energy eigenvalues. 
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FIG. 1: Andreev levels for a s± superconductor S/I/S junc- 
tion with r — 2 and increasing values of a. (a): Z = 0.7. 
Black line, a = 0.01, red line, a = 0.5, green line, a = 0.7. 
(b): Z — 2. Black line, a — 0.01, red line, a — 0.4, green line, 
a — 0.7. The black line curves correspond to a nearly perfect 
(a = 0) "conventional" s-wave junction. The blue lines are 
the Andreev level for a — a c — 0.7071. 

In the nearly insulating limit (Z — >• 00), the system 
decouples and we obtain information on the two sepa- 
rate electrodes. More precisely, in this limit, the spectral 
equation (|6|) describes surface bound states of energy 

/l -r 2 a 4 

already discussed in ref. [8| for a junction N/I/s ± . 



For < a < a c , the two emerging levels E 2 = 
±62(^5 ol) and the levels E\ — d=ei(y?, a) have, in gen- 
eral, opposite dispersions, i.e. dE 2 /d(p > 0, dEi/dcp < 
as can be seen in Fig. [Ha) and (b). The sign difference 
of the dispersions is the key point in determining the 
Josephson current-phase relation, as it will be discussed 
below. 

The Josephson current Id carried by the discrete An- 
dreev levels Ek through the contact can be found from 
the free energy, according to the following relation [28|, [3l| 

^W^tanh^J) (10) 

where f(Ek) is the Fermi distribution function and fc 
labels the Andreev level with energy Ek = ±efe(<^, a). 
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FIG. 2: current-phase relation for a s± superconductor S/I/S 
junction with r = 2, T/T c = 0.01 and increasing values of a. 
(a)Z = 0.7, (black line, a = 0.01, red line, a = 0.5, green line, 
a = 0.7). (b) Z = 2. (black line, a = 0.01, red line, a = 0.4, 
green line, a = 0.7) 

Fig. [2](a) and (b) show the current-phase relations 
Id(<p,ct), at low temperature, corresponding to the dis- 
crete spectrum of Andreev levels represented in Fig. [1] 
(a) and (b) respectively and calculated through Eq. ([TP]) . 
The current in this figure is normalized with respect to 
Io = eAi(0)/fi, which is the zero temperature maxi- 
mum Josephson current through a one dimensional s- 
wave junction in the clean limit. The three curves corre- 
spond to values of the band ratio approaching the crit- 
ical value a c = 0.7071. The current-phase relation rep- 
resented by the green lines (a = 0.7) in both Figs [2] 
(a) and (b) shows a clear tt phase-shift (the maximum 
Josephson current for < cp < tt has a negative value). 
For the considered a value this tt state persists in the 
whole temperature range (see Fig. [3] (a) and (b)). The 
mechanism of formation of this state is the following. 
As the band ratio increases an upper +c 2 ((p, a) and a 
lower —c 2 ((p, a) extra Andreev bands gradually emerge. 
These two bands have different character compared to 
the low energy bands d=ei(y?, a): they transport super- 
currents in the opposite directions. As the temperature 
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decreases, only low energy level are populated while those 
at higher energy are empty. In the competition between 
these opposite carrying current energy levels ±62 (y?, a) 
and ±ei(<^, a), which coexist for any value of the band 
ratio, it is the temperature that determines the direction 
of the total current and the possible existence of a tt — » 
crossover in the considered two gap superconductor one- 
dimensional junction. 

The situation is similar, although not identical, to that 
found in a two-dimensional d-wave tt junction. In this 
case [24], two kinds of bands, conventional and midgap, 
alternate, without coexisting and compete each other in 
determining the supercurrent, depending on the angle of 
incidence of the Andreev quasiclassical trajectories with 
the interface. Therefore it is the two-dimensionality here 
that plays the key role. 

In Fig. (j3]), (a) and (b), we report, for increasing val- 
ues of a , the dependencies of the normalized Josephson 
critical current I c /Iq (I c — m ax^/^((^)) from the reduced 
temperature T/T c , derived by the discrete Andreev spec- 
tra for a nearly clean junction [Z = 0.001) and a low 
transparency junction (Z = 3), respectively. 

The negative sign of I c (blue lines, Fig. [3] a) and b)) 
indicates that the junction free energy, i.e. the quantity 
F((p) = ^o/2tt d(pl(ip), has a minimum at <p = 7r such 
that the ground state of the junction is tt shifted. As 
discussed above, the tt — > crossover occurring with de- 
creasing temperature for a = 0.65 in figure [3] (a) and for 
a = 0.6 (green line) in figure [3] (b)(see inset) is due to 
the competition between the contribution to the current 
from the high energy band E2 and that low energy band 
El 

So far, by using Eq. ([T0]h we have included the contri- 
bution to the Josephson current from the discrete energy 
spectrum only. However the continuous-spectrum states 
make their own contribution to the current which has to 
be accounted for properly. Following the approach de- 
veloped by Furusaki Tsukada [32[ the total Josephson 
current, including contributions from both the Andreev 
bound states and the continous spectrum is given by 

eAi/c B T 1 r / • \ / • \i /1 1 \ 

a{-(p,icj n ,a)\ (11) 



I 



n 



and introduced 



Here we have defined ^1 = y^n 
the Matsubara frequencies uo n = 7r/c#T(2n +1), with n 
ranging from —00 to +00. a(ip, iuj ni a) is a scattering am- 
plitude coefficient for the process in which an electron- 
like quasiparticle traveling from the left of the junction 
is reflected back as a hole-like quasiparticle in the pres- 
ence of two bands. This coefficient is derived by solving 
the BdG equations under the assumption E > Ai and 
dropping the requirement of exponentially decay of the 
solutions for \x\ — > 00. The details of this procedure will 
be given elsewhere [33|. The results of the calculated 
total critical current I c as a function of the reduced tem- 
perature, are shown in Fig. 2J (a) and (b), for different 
values of the band ratio a and already considered in Fig. 
® 
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FIG. 3: Critical current as a function of temperature for a 
s± superconductor S/I/S junction with r = 2, T/T c = 0.01 
(a)Z = 0.001; black line, a = 0.001, red line a = 0.5, green 
line, a = 0.65, blue line, a = 0.7. (b) Z = 3; black line, 
a = 0.001, red line, a = 0.5, green line, a = 0.6, blue line, 
a = 0.7. The negative sign of I c indicates that the junction 
minimum is at ip = 7r (7r-j unction). The inset in panel (b) 
shows, on a larger scale, the transition ty — » for a — 0.6 
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FIG. 4: Critical current as a function of temperature. The 
contribution of the continuous energy spectrum states has 
been added. (a)Z = 0.001;(b)Z = 3; in both figures, (a) 
and (b), the black, red, green, brown, blue lines correspond 
to the curves with a = 0.001, 0.5, 0.6, 0.65, 0.7, respectively. 
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The general features of the curves I c vs T, shown 
in Fig. (J2J) qualitatively reproduce those of Fig. (|3]). 
Most notably they confirm the ir — junction crossover. 
However some noticeable differences may be pointed out. 
For instance, the tt — step-like crossover, calculated for 
a = 0.65 and shown in Fig. (j4]) (a), is smoothed out when 
considering the continuous spectrum contribution, as can 
be seen by comparing with the corresponding curve in 
Fig. ©, (b). 

In conclusion the model investigated in this paper for 
a multiband superconductor symmetric s ± Is ± junction, 
predicts a number of non trivial details related to the 
Josephson effect in these systems. These results are of 
relevant interest for the case of all-pnictide Josephson 
micro-junctions. The spectral equation provides the An- 
dreev bound levels as a function of the band ratio param- 
eter a. The main effect of the presence of the second band 
is the building up of two extra Andreev levels which drive 
Cooper pairs in a direction opposite to that observed in 
the presence of a single band. 

The phase-current relations predicted on the basis of 



the Andreev levels (without consideration for the contin- 
uous energy spectrum) shows the formation of a note- 
worthy 7r-state in the junction coupling. For selected 
values of the band ratio, a tt —> crossover may occur 
as the temperature decreases. In this case the 7r-state 
is observed at high temperature whereas the 0-state is 
observed below the crossover temperature. In this re- 
gion the junction recovers the behavior of a conventional 
"0" junction. These results are confirmed by means of 
a more exhaustive evaluation of the Josephson current 
englobing the contribution of the continuous spectrum 
energy states. 
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